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Abstract 

We present the preliminary result of a search for the decay of K*^p~. The data were 

recorded with the BABAR detector at the PEP-II collider and correspond to 123 million BB pairs 
produced in the e~^e~ annihilation through the T(45) resonance. We obtain an upper limit on 
the branching fraction for this decay of B{B^ —>■ K*^p^) < 24 x 10^^ (90% C.L.). Ah results are 
preliminary. 
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1 INTRODUCTION 



The study of the B meson decay into charmless hadronic final states plays an important role in the 
understanding of the origin of CP violation. The B decays to two vector particles are of special 
interest because their angular distributions reflect both strong- and weak-interaction dynamics pQ. 
The decay B — > K*p is dominated by 6 — > s penguin contribution, and the tree-level contribution is 
CKM suppressed in the Standard Model (SM) (see Figure^. The angular correlation measurement 
is particularly sensitive to phenomena beyond the SM, potentially present at either loop- or tree- 
level 0. 




Figure 1: Gluonic penguin and tree diagrams contributing to the process B K*^ p~ . The 
penguin contribution coming from the diagrams with t and c quarks in the loop dominates since 
contributions from process with a u quark is suppressed. The tree is a CKM-suppressed diagram. 

For measurements with limited statistics, we integrate over the azimuthal angle (j) between the 
two decay planes shown in Figure El The differential decay rate is 

Tdcoselldcosep- = I (^^ cos'(^K*^cos%- + \{1 - fL)s^n'e^,,s^n%-^ (1) 

where is the longitudinal polarization fraction component /l = F^/F [Hll^. The angles 0^*+ ^^'^ 
9p- are the helicity angles of K*'^ and p^, which are defined between the charged K['k) direction 
and the direction opposite the B in K*^[p^) rest frame as shown in Figure [21 




Figure 2: Definition of angles 6k*+, &p-, and 0, for the decay B^ K*^p . The K^-k^ (vr tt") 
final states are shown in the K*~^ ip~) rest frame. 

In the SM, and assuming naive factorization, the polarization is expected to be proportional 
to (1 — 4 X /m|) > 90% 12], which has been verified experimentally in both B^ P^P^ and 
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— > p~^p decays jSlinilZl- However, this prediction does not agree with measurements in pure 
pengum B decays such as B+ (j)K*+ p and B° (j)K*^ |H] as shown in Table □ Since the pure 
penguin B decay processes are sensitive to new physics, it is very interesting to look at additonal 
pure penguin or penguin dominated B decays |2]. 

Table 1: Previous measurements for rates and /l for B ^ V1V2 modes. The first error is statistical 
and the second error is systematic. 



Mode 


fi(io-**) 


h 


B+ ^ 




22.5l^:^±5.8 


0.97+^:^^ ± 0.04 


5°^ 




30±4±5 


0.99 ± 0.03 ± 0.04 


B+ ^ 
5°^ 


p^K*+ 


10.6^1^ ± 2.4 
9.2 ± 0.9 ± 0.5 


OMtoit ± 0.04 
0.52 ±0.05 ±0.02 


B+ ^ 


(j)K*+ 


12.7i|^ ± 1.1 


0.46 ±0.12 ±0.03 



In this paper we report a search for the B" K*^p decay based on a sample of 123 million 
BB pairs and set a limit on the branching fraction. 

2 THE BABAR DETECTOR AND DATASET 

The results presented in this paper are based on data collected in 1999-2003 with the BABAR detec- 
tor |HI at the PEP-II asymmetric e~^e~ collider located at the Stanford Linear Accelerator Center 
(SLAC). An integrated luminosity of 113 fb~^, corresponding to 123 million BB pairs, was recorded 
at the T{4:S) resonance ( "on-resonance" ) with the center-of-mass(CM) energy ^/s = 10.58 GeV. An 
additional 12 fb~^ were taken about 40 MeV below this energy ("off-resonance") for the study of 
continuum backgrounds in which a light or charm quark pair is produced instead of an T(4S'). 

The asymmetric beam configuration in the laboratory frame provides a boost of /?7 = 0.56 to 
the T(4S). Charged particles are detected and their momenta measured by the combination of 
a silicon vertex tracker (SVT), consisting of five layers of double-sided silicon strip sensors, and 
a 40-layer central drift chamber, both operating in the 1.5-T magnetic field of a solenoid. The 
tracking system covers 92% of the solid angle in the CM frame. 

Charged-particle identification (PID) is provided by the average energy loss (dE/dx) in the 
tracking devices and by an internally refiecting ring- imaging Cherenkov detector (DIRC) covering 
the central region. A K/tt separation of better than four standard deviations (cr) is achieved for 
momenta below 3 GeV/c, decreasing to 2.5 a at the highest momenta reached by the B decay final 
states. Photons and electrons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). The 
EMC provides good energy and angular resolutions for detection of photons in the range from 30 
MeV to 4 GeV. The energy and angular resolutions are 3% and 4 mrad respectively, for a 1 GeV 
photon. The flux return for the solenoid is composed of multiple layers of iron and resistive plate 
chambers for the identification of muons and long-lived neutral hadrons. 
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3 ANALYSIS METHOD 



Monte Carlo (MC) simulations ^Hl of the signal decay modes, continuum and BB backgrounds are 
used to establish the event selection criteria. We reconstruct B^ K*^p~ candidates from the 
decay products of the K*~^ K~^ir^ and p~ 7r~7r*^ (the charge conjugate states are implied in 
this paper). Charged-track candidates are required to originate from the interaction point: distance 
of closest approach to the interaction point less than 10 cm along the beams direction and less than 
1.5 cm in the plane transverse to the beams direction. We require that the track from the p~ decay 
has particle identification information inconsistent with the electron, kaon, or proton hypotheses; 
and the track from the K*~^ decay should agree with the kaon hypothesis. We reconstruct 
mesons from pairs of photons, where each photon must have an energy greater than 50MeV in 
the laboratory frame and must exhibit a lateral profile of energy deposition in the electromagnetic 
calorimeter consistent with an electromagnetic shower 0. The vr'^ candidates must have a mass 
that satisfies 0.11 < m{-f-f) < 0.16GeV/c^ The mass of the reconstructed p and K*^ candidates 
must satisfy 0.396 < m(7r-7r0) < 1.146 GeV/c^ and 0.767 < m{K+7r^) < 1.017 GeV/c^, respectively. 
Combinatorial backgrounds dominate near 100,3(0^)1 = 1 (V denotes K*~^ or p~). Backgrounds 
from B decays, like B^ — > K*^{892)7r~, with an additional low energy vr^ from the rest of the 
event (ROE), tend to concentrate at negative values of cos{6v)- The K*~^ K~^7r^ and p~ 
TT'ir^ helicity angles are restricted to the region —0.8 < cos{9y) < 0.98 to suppress combinatorial 
background and reduce acceptance uncertainties due to low-momentum pion reconstruction. 

The i?-meson candidates are identified from two nearly independent kinematic observables [S], 
the beam energy-substituted mass mES = [(■s/2 + Pi • PB)'^/Ef — p'^]^^'^ and the energy difference 
AE = (EiEs — Pi ■ Pb — s/2)/^/s, where {Ei,pi) is the four-momentum of the e~^e~ initial 
state, and {Eb,Pb) is the four- momentum of the reconstructed B candidate, all defined in the 
laboratory frame. For signal events, the mES distribution peaks at the B mass and the AE 
distribution peaks near zero. We accept candidates that satisfy 5.21 < niES < 5.29 GeV/c^ and 
—0.12 < AE < 0.15 GeV. The asymmetric AE window suppresses background from higher- 
multiplicity B decays. 

In this analysis, B^ decays to charm modes, such as B^ — > D'^ir^ with — > K~^7r~7r^ have the 
same final state as signal. If the tracks from these B decays are used to reconstruct the K*^ and 
p~ mesons, these events will have peaking AE and mEs distributions under the signal region. We 
apply the requirements |m(K+7r") - m{D^)\ > 0.02 GeV/c^ and |m(/C+7r"7r°) - m{D^)\ > 0.04 
GeV/c^ to suppress these peaking backgrounds, where m{D^) is the nominal mass of meson |12) . 
After veto 95% (99%) of longitudinal(transverse) signal events are retained. 

Signal candidates may pass the selection even if one or more of the tracks or vr^s assigned to the 
K*^p~ state actually comes from the other B in the event. These self-cross-feed (SCF) candidates 
comprise 37% (21%) of the accepted signal for longitudinal (transverse) signal and are included as 
signal in the fit. 

Continuum e~^e~ qq {q = u,d,s,c) events are the dominant background. To discriminate 
signal from continuum we use a neural network (NN) to combine six variables: the Fisher of the 
Legendre monomials jll| : the sum of transverse momenta in the ROE relative to the z axis; the 
cosine of the angle between the direction of the B and the collision axis (z) in the CM frame; the 
cosine of the angle between the B thrust axis and the z axis; the cosine of the angle between the 
S-thrust axis and the thrust axis of the ROE; the decay angle of one of the vr'^s (defined in the 
same way as the K* / p decay angle, Oy)-, randomly selected. The final sample of signal candidates 
is selected with a cut on the A^A^ output that retains ~ 93% (54%) of the signal (continuum). 
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When multiple B candidates can be formed, we select the one that minimizes the sum of the 
of the reconstructed vr^ masses from the nominal mass. For those with the same lowest x^; 
we keep the first one. 

The efficiency of the selection is 6.8% (13.9%) for longitudinally (transversely) polarized signal 
as determined with MC simulations. After the full selections, we obtain 14251 events in the data 
sample, which are dominated by combinatoric backgrounds: roughly 92% from qq and 7.7% from 
BB. 

We use MC-simulated events to study the cross- feed from other B decays. The charmless modes 
are grouped into thirteen classes with similar kinematic and topological properties. Two additional 
classes account for the neutral and charged 6 — > c decays. For each of the background classes, a 
component is introduced into the likelihood below, with a fixed number of events. In the selected 
K*~^p~ sample we expect 56 ± 27 charmless background events and 1005 b ^ c events. 

We use an unbinned, extended maximum-likelihood (ML) fit to extract the signal yield. The 
likelihood for each B candidate i is defined as 

TV . . 

= e-^' J] Ar-9p-9 + ^gcipm + ^ ^.^B I (2) 
i=i ^ j=i ^ 

where N' is the sum of the signal and continuum yields and the fixed S-background yields, A^**^ is 
the number of signal events of type K*^p^ in the entire sample. iV?^ is the number of continuum 
background events and is floated in the fit. The numbers of events rij in the B background 
category j are all fixed to their MC expectations. The probability density function (PDF) V 
is the product of the PDFs of seven discriminating variables. The signal PDF is thus given by 
psig ^ v[rnEs) ■ V{AE) ■ V{NN) ■ V{cos{9k*+)) ■ V{Mvi) ■ V{cos{dp-)) ■ V{Mv2), where Myi 
{Mv2) is the mass of K*^ {p~) meson, and V{cos{9x*+)), 'P{cos{6 p-)) is the signal helicity PDF 
which is expressed as a function of the longitudinal polarization (see Eq. The ideal angular 
distribution is multiplied by the detector acceptance function. We obtain the acceptance function 
from a simultaneous fit to a sample of MC events with transverse and longitudinal polarization. 
The PDF of the continuum contribution is denoted V^^. The corresponds to the PDF of the B- 
background category j. The signal events are decomposed into two parts with distinct distributions: 
signal events that are correctly reconstructed and those mis-reconstructed, namely, SCF events. 
The SCF fractions for longitudinal and transverse signal are estimated by MC simulation. The 
niES-: NN, cos{9x*+), Myi, cos{6p-) and Mv2 PDFs for signal and B background are taken 
from the simulation. The continuum-background mESi A£^, NN, cos{9j^*+) and cos{6p-) PDF 
parameters are fioated in the fit to data. The distributions of the continuum as a function of Myi 
and My2 are described by a non-parametric PDF '13' derived from mEs and IS.E data sidebands. 
A total of 12 parameters, including signal yield and continuum background yield, are varied in the 
fit. 

4 SYSTEMATIC UNCERTAINTIES 

The contributions to the systematic error on the signal parameters are summarized in Table |21 
The uncertainties due to the signal model are obtained from varying the signal PDF parameters, 
which are fixed in the fit, within their estimated errors and assign the effects on the signal yield as 
systematic error. We perform fits on large MC samples with the measured proportions of K*^p~ 
signal, continuum and B backgrounds. The bias observed in these tests is due to imperfections of 
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the PDF model: e.g., unaccounted correlations between the discriminating variables of the signal 
and B-background PDFs. The bias is assigned as a systematic uncertainty of the fit procedure. 
The expected event yields from the B background modes are varied according to the uncertainties 
in the measured or estimated branching fractions. 

Table 2: Summary of the systematic uncertainties in the measurement of the — > K*~^p~ branch- 
ing fraction. 



Source 


Uncertainty 


Fit uncertainties (in Events) 


Signal model 


+2.3 
-1.9 


Fit procedure bias 


4.2 


B backgixMuidti 


1.7 


Total fit error 


5.1 


Multiplicative [%] 


Track finding 


2.4 


Neutral correction 


10.3 


Number of BB pairs 


1.1 


Particle ID 


1.1 


Total multiplicative uncertainties 


10.7% 


Non-resonant charmless background 


+U.U w 
-39.0 



In this analysis, we do not include a fit component for other B decays with the same final- 
state particles selected within the K* or p resonance mass window, such as the non-resonant 
decays B^ — ^ K~^'!T~ir'^7r^^ , B^ K'^tt^ and B^ K*^n^7r^\ The selection requirements alone 
suppress the B^ K^tt^tt^tt^ [B^ p~ K~^tt^ and B^ K*^^^^^) efficiency by two (one) orders 
of magnitude relative to B^ — > K*~^p~ . The contribution of these decays to the fit results is also 
significantly suppressed by the mass and helicity-angle information in the fit; they are examined in 
the context of mass and helicity-angle distributions, as discussed below. 

To check the sensitivity of our results to the presence of non-resonant B^ K~^TT~-K^7r^ , B^ — > 
p~ K^TT^ and B^ —>■ K*~^Tr~Tr^ decays, we explicitly include a fit component for them, assuming 
a phase-space decay model. The associated systematic error is estimated by the difference in the 
data fit result when the yields of these background modes are floated or fixed to zero. We obtain 
an asymmetric error of -22 events (-39%) on the signal yield, which is systematically overestimated 
when these non-resonant background modes are not modeled in the ML fit. This systematic error is 
preliminary estimation, and is presented separately from the other systematics, with the label "non- 
resonant". Interference effects between the resonant and non-resonant components are ignored in 
this fit. 

The systematic uncertainties in the efficiency arc due to track finding (2.4% for two tracks), 
particle identification (1.1% for two tracks), and vr*^ reconstruction (10.3% for two tt^s). Smaller 
systematic uncertainties arise from event-selection criteria, MC statistics, and the number of B 
mesons in the sample. 
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5 PHYSICS RESULTS 



From the ML fit, we find the signal yield Ngig = 58 it 19(siat). The results are summarized in 
Table |21 It is checked that we obtain the most conservative estimate for the upper limit on the 
branching fraction when 70% longitudinal polarized signal is used in the fit. Using the above results, 
together with the selection efficiency, the branching fractions of B{K*^ — > K^t:^), B{p~ — > 7r~7r'^), 
B{'K^ — > 77), we obtain a central value of the branching fraction, B{B^ K*^p^) = [16.3 it 
5.4(stat) lb 2.3(syst)lg'3(non — resonant)] x 10~^. The impact of the uncertainties on B{K*^ — > 
K^TT^), B{p~ — > 'K~-K^), B{tt^ 77) is negligible compared to the other systematic errors. 

In Table |31 the statistical error on the branching fraction results from the statistical errors on 
the signal yield. The systematic uncertainty on the branching fraction results from the propagation 
of the systematic uncertainties on signal yield and effective selection efficiency. Finally, the non- 
resonant systematics on signal yield is propagated to the branching fraction. 

Figure 01 shows the result of the fit projected onto the tues, ^E, Myi and Mv2 observables. 
The histograms show the data after a cut on the quantity Vsig/ i'Psig + 'Pcont) has been applied, 
where Vsig and Vcont are the probabilities for a given event to be signal and continuum background, 
respectively, and are evaluated using all observables except for the one that is being plotted. 

Because of the limited statistical significance of the observed signal, we choose to quote as our 
preliminary result an upper limit on the branching fraction. Taking systematic uncertainties into 
account, an upper limit of 24 x 10^^ at 90% confidence level (C.L.) is set for the — > K*^p~ 
branching fraction. 

Table 3: Summary of the fit results: Signal yield (Ngig), effective selection efficiency (e), branching 
fraction (B), upper limit at 90% confidence level and significance of the measurement, expressed 
as number of standard deviation (crs). The first error corresponds to the statistical uncertainty 
and the second one to the systematic uncertainty, and the third is the systematic uncertainty from 
non-resonant contributions. 



Quantity 


Measured Value 


Nsig 


58 ± l%stat) 


e{%) 


8.9 ± 1.0 


S(xl0-6) 


16.3 ± 5.4(stat) ± 2.3(syst) tgislnon - resonant) 


C/.L.(xlO-6) 


24 (22 statistical only) 


Significance (a) 


3.2 (3.7 statistical only) 



6 SUMMARY 

We have searched for the decay K*^p^ using a maximum likelihood technique in a data 

sample equivalent 113 fb^^ of integrated luminosity. From a fitted signal yield of 58 it 19(stat), we 
obtain a branching fraction is B{B^ K*~^ p~) = [16.3ib5.4(stat)ib2.3(syst)]l^g3(non — resonant)] x 
10~^. We get a preliminary upper limit on the B^ K*^p~ branching fraction at 90% C.L. is: 

B{B^ K*+p-) < 24 X 10"^ 

All results are preliminary. 
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Figure 3: The (a) mEs, (b) AE, (c) K* mass, (d) p mass distribution for signal enriched samples 
of the data. The dashed line is the projection of the continuum background, the dotted line is the 
projection of the sum of backgrounds and the solid line is the projection of the fit result. 
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